PD-A184  888  5TUDV  OF  BOUNDARV  STRUCTURES (U>  WASHINGTON  UNIV  SEATTLE  1/1 
"  DEPT  OF  MATERIALS  SCIENCE  AND  ENGINEERING  R  KIKUCHI 

31  JUL  87  ARO-19984  3-PH  DAAG29-85-K-8166 
UNCLASSIFIED  F/G  11/6  NL 


fiflo  IWoH  S  'Pti 


STUDY  OF  BOUNDARY  STRUCTURES 


FINAL  REPORT 


RYOICHI  KIKUCHI 


JULY  31,  1987 

U.S.  ARMY  RESEARCH  OFFICE 

DAAG29-85-K-0166  (U.W.) 
DAA629-83-C-0011  (HUGHES) 


DEPARTMENT  OF  MATERIALS  SCIENCE  AND  ENGINEERING 
UNIVERSITY  OF  WASHINGTON 
SEATTLE,  WASHINGTON  98195 


OTIC 

£p|ELECTE 
V\  SEP  1  7  1987 

APPROVED  FOR  PUBLIC  RELEASE;  ^  CW 
DISTRIBUTION  UNLIMITED. 


87  9 


<)  095 


THE  VIEW,  OPINIONS,  AND/OR  FINDINGS  CONTAINED  IN  THIS  REPORT  ARE  THOSE 
OF  THE  AUTHOR(S)  AND  SHOULD  NOT  BE  CONSTRUED  AS  AN  OFFICIAL  DEPARTMENT 
OF  ThE  ARMY  POSITION,  POLICY,  OR  DECISION,  UNLESS  SO  DESIGNATED  BY 
OTHER  DOCUMENTATION. 


la  REPORT  SECURITY  CLASSIFICATION 


REPORT  DOCUMENTATION  PAGE 

RESTRICTIVE  MARKINGS 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 

2b  OECLASSlFlCATlON  /  DOWNGRADING  SCHEDULE 

4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 
Requisition/Purchase  Request 
P-23188-PH  U.W. 


|6b  OFFICE  SYMBOL 
(If  applicable) 


6a.  NAME  OF  PERFORMING  ORGANIZATION 
University  of  Washington 
Dept,  of  Mat.  Sci.  &  Engr. 

6c.  ADDRESS  (O' ty,  State,  and  ZIP  Code) 
Mail  Stop  FB-10 
Seattle,  WA  98195 


8a.  NAME  OF  FUNDING /SPONSORING  8b.  OFFICE  SYMBOL 

ORGANIZATION  (If  applicable) 

U.  S.  Army  Research  Office 

8c,  ADDRESS  (City,  State,  and  ZIP  Code) 

P.  0.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 

1 1  TITLE  (Include  Security  Classification) 

Study  of  boundary  structures  (Unclassified) 

12  PERSONAL  AUTHOR(S) 

Ryoichi  Kikuchi _ 

1 3a  TYPE  OF  REPORT  13b  TIME  COVERED 

Final  Report  from  4/1/83  to 7/31/87 

16  supplementary  notation 


3  distribution /availability  of  report 

Approved  for  public  release; 
distribution  unlimited. 

5  MONlTOB'NG  ORGANIZE  T,^N  REPORT  NUMBER(S) 

ARO  19904. 3-PH 
ARO  23188. 1-PH 

7a.  NAME  OF  MONITORING  ORGANIZATION 

U.  S.  Army  Research  Office 
7b.  ADDRESS  (City,  State,  and  ZIP  Code) 

P.  0.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 

9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

DAAG29-83-C-001 1  DAAG29-85-K-0166 

10  SOURCE  OF  FUNDING  NUMBERS  ” 


PROGRAM 
ELEMENT  NO 


PROJECT 

TASK 

NO 

NO 

WORK  UNIT 
ACCESSION  NO 


14  DA"-  nt  ocport  (Year,  Month,  Day)  15  PAGE  COUNT 

July  31,  1987 


The  view,  opinions  and/or  findings  contained  in  this  report  are  those 


1  7 

COSATI  COOES  i 

FiELO 

GROUP 

SUB-GROUP 

18.  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 


'9  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

The  project  extends  W  previous  study  on  structures  of  a  2-D  grain  boundary  (g.b.). 
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The  quadruple  square  (QS  treatment.  A  square  lattice  is  chosen  as  the  basic  lattice 
(called  the  DSC  lattice)  of  the  analysis.  The  first,  second,  and  third  nearest  neighbor 
pairs  are  excluded  due  to  a  large  repulsion.  The  fourth  and  fifth  neighbors  contribute 
interaction  potentials  and  ,  a^d  so  on.  In  the  DSC  lattice  the  CVM  is  formulated 

using  a  cluster  made  of  'nine  lattic^1 points  in  the  form  of  a  quadruple  square.  The  QS  form¬ 
ulation  leads  to  numerical  results  that  are  reliable  compared  with  the  MC  method. 

The  pair  treatment.  In  low  temperatures,  atoms  hardly  deviate  from  the  skeletal  stable 
crystalline  structure^  and  hence  we  can  ignore  DSC  lattice  points  not  on  the  stable  crystal 
which  is  a  square  lattice.  Then  we  can  use  the  pair  cluster  as  the  basic  cluster.'  The  pair 
method  is  found  to  give  reliable  numerical  results  practically  the  same  as  those  of  QS  for 
low  temperatures  including  the  entire  ordered  phase  in  the  work  D.  This  result  is  signifi¬ 
cant  because  when  the  pair  method  is  used,  the  two  sides  of  the  g.b.  do  not  need  to  be 
oriented  symmetrically  as  were  done  in  the  present  report.  Applications  of  this  point  of 
view  for  g.b.  of  different  angle  combinations  and  the  3-D  cases  are  now  planned. 

Hybrid  cluster.  In  the  work  D  for  the  ordered  system  g.b.,  it  is  found  necessary  to 
supplement  the  pair  cluster  by  a  triangle  cluster  near  the  coincident  lattice  site.  Such 
a  hybrid  choice  of  clusters  was  tried  for  the  first  time,  but  it  makes  sense.  We  will  need 
more  hybrid  choices  near  the  g.b.  center  in  order  to  take  into  account  the  frustration 
effect  of  atoms. 


The  low  temperature  "melting"  of  g.b.  In  all  of  A,  B,  and  D  cases,  we  found  that  the 
low  T  behavior  and  the  high  T  behavior  of  g.b.  can  be  clearly  distinguished.  In  low  T,  it 
is  crystal-like.  As  T  increases,  the  g.b.  opens  up  as  -log  (T  -  T),  with  T  being  the 
bulk  melting  temperature.  The  transition  region  between  the  two  is  about  one-half  of  T^. 

Diffusion  along  the  g.b.  Diffusion  along  the  g.b.  is  orders  of  magnitude  larger  near 
the  g.b.  than  the  bulk  diffusion.  The  diffusion  coefficient  is  proportional  to  the  square 
of  the  local  vacancy  density.  This  holds  when  the  impurities  are  adsorbed  to  the  g.b.  and 
also  when  they  are  repeated. 
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1.  Introduction  and  Summary 


In  any  type  of  crystalline  materials,  he  it  metallic,  ceramic,  dielectric,  or  semiconducting,  grain  bound¬ 
aries  (g.b.)  play  the  key  roll  in  affecting  mechanical  properties  (strength  and  toughness).  However,  the 
atomic  structure  of  g.b. has  not  been  known,  except  the  T~0  behavior,  mainly  because  of  the  lack  of 
techniques  which  can  analyze  the  structure  The  present  project  studies  the  structure  anil  properties  of 
g.b.  in  a  two-dimensional  system  for  finite  temperatures  using  statistical  mechanical  techniques  developed 
by  the  author. 


In  one  of  the  tasks  performed  under  the  contract  preceding  the  present  one,  a  2-D  grain  boundary 
(g.b.)  was  studied  theoretically.1  The  free  energy  /'of  the  entire  system  (including  the  g.b.)  was  written 
in  terms  of  the  variables  (called  the  state  variables)  which  represent  probabilities  of  local  configurations. 
The  special  feature  of  this  work  was  that  the  cluster  variation  method  (CVVIp  was  used  in  formulating 
the  entropy  in  terms  of  the  basic  variables.  By  minimizing  /',  the  equilibrium  state  of  the  g.b.  was  derived, 
including  the  density  profile  and  the  excess  free  energy  value.  Since  this  work  forms  the  basis  of  the  tasks 
done  under  the  contract  being  reported  here,  we  attach  Ref.  I  as  Appendix  A. 

One  of  the  main  conclusions  of  Ref.  1  was  that  the  g.b.  starts  to  melt  at  around  one-half  of  the 
mcltine  temperature  /„,,  and  the  g  b.  broadens  near  / as  -  !og(  /„,  -  /').  This  conclusion  attracted  the 
interest  of  workers  and  several  papers  have  been  published  on  the  subject  reporting  the  Monte  Carlo 
simulations.  The  low  temperature  g.b.  melting  was  confirmed  by  the  simulations,  the  latest  paper  of  which 
may  be  the  one  due  to  Yip  et  al.1 


Since  the  method  used  in  Ref.  1  is  versatile  in  studying  properties  of  grain  boundaries,  we  extended 
it  to  several  cases  during  the  present  contract.  The  results  are  in  Appendices  B.  C,  and  D,  and  they  are 
summarized  in  subsequent  sections. 


Appendix  B  is  an  extension  of  Ref.  I  and  includes  impurities.  Appendix  C  is  based  on  B  and 
calculates  atomic  diffusion  along  the  g.b.  This  paper  is  significant  in  that  the  technique  of  irreversible 


statistical  mechanics  (the  pat li  probability  method’1  PPM)  is  used  in  the  g.b.  studies  It  is  shown  that 
diffusion  is  orders  of  magnitude  larjicr  near  the  g.b  compared  with  the  bulk  diffusion.  One  significant 
result  is  that  for  either  case  when  the  impurities  are  adsorbed  to  nr  rejected  from  the  g.b.,  the  diffusion 
coefficient  becomes  large  near  the  g.b.  in  proportion  to  the  square  of  the  local  vacancy  concentrations. 


Appendix  D  works  with  the  g.b.  be* ween  two  ordered  phases.  [  he  symmetric  boundary  and  asym¬ 
metric  boundary  are  distinguished.  I  lie  profile  of  the  long-ianue  order  across  the  itb  is  calculated. 


In  the  boundary  problems  the  number  of  independent  variables  is  usually  large,  easily  exceeding  |('()0. 
and  a  minimum  of  the  free  energy  (actually  the  grand  potential  if)  is  to  be  found  with  respect  to  these 
gambles  1 1,  •  linear  iteration  technique  called  the  NIM  works  amazingly  well  and  Li  always  decreases 
as  the  iteration  proceeds,  however  small  the  decrement  may  be. 


I  he  importance  of  the  gram  boundaries  in  influencing  mechanical  properties  of  materials  comes  mainly 
from  the  concentration  of  vacancies  at  the  g  b.  Statistical  mechanics  can  calculate  the  amount  of  vacancies 
at  the  g  h  as  we  have  done  in  this  project  and  reported  in  the  appendices.  The  concentration  ol  vacancies, 
however,  is  only  the  beginning  of  the  g  b  studies  of  practical  importance.  It  is  needed  to  proceed  to  study 
how  the  i  .n  ancies  influence  kinetic  be.h.n  ior  of  the  materials. 


I  he  diffusion  study  in  Appendix  ('  is  one  step  ahead  o.  the  prevailing  effort  of  statistical  mechanics, 
but  'till  is  far  short  of  the  really  important  aspec  t  of  the  kinetic  study,  namely  the  plastic  deformation  and 
fracture  It  is  planned  that  we  w ill  mm  c  on  to  the  study  of  grain  boundary  contribution  to  these  phenomena. 


I  he  CVM  and  PPM  are  \  cry  yvcll  suited  in  study  mg  the  structure  ol  boundaries  analy  tically.  Although 
the  yyay  may  still  be  long  to  achieve  the  goal  of  understanding  mechanical  properties  ol  materials  from 
microscopic  behavior,  we  will  continue  our  effort  towards  the  goal  by  improving  these  methods. 


Introduction  and  Summary 


2.  Grain  boundaries  with  impurities  in  a 
two-dimensional  lattice  gas  model  (attached  as 

Appendix  B) 


I'hc  model  of  Ref.  1  was  extended  to  include  impurities.  The  work  was  accepted  for  publication  in 
Physical  Review  B,  and  is  attached  as  Appendix  B  in  this  report.  I'hc  g.b.  is  formed  on  the  basic  DSC 
lattice.  In  the  CVM  analysis  we  use  a  quadruple  square  (QS)  of  nine  DSC'  points  as  the  basic  cluster. 


When  the  matrix-impurity  interaction  is  repulsive  so  that  the  impurities  arc  rejected  from  the  bulk 
matrix,  they  arc  adsorbed  near  the  grain  boundary'.  This  is  because  vacancies  are  concentrated  around  the 
g  b.  When  the  matrix  and  impurity  interaction  is  attractive,  impurities  stay  away  from  the  g.b.  The 
amount  of  adsorption  is  calculated  as  a  function  of  temperature. 


As  was  observed  in  Ref.  1,  the  low  temperature  structure  and  the  high  temperature  structure  of  the 
g.b.  can  be  clearly  distinguished.  In  the  high  temperature  region,  the  g.b.  widens  as  in  Ref  1,  depending 
on  I  as  -  log( /„,  -  7).  This  behavior  is  derived  based  on  the  calculations  of  both  the  entropy  and  the 
adsorption. 


2  ( irain  boundaries  with  impurities  in  a  two-dimensional  lattice  gas  model  (attached  as  Ap¬ 
pendix  H) 


3.  Theory  of  diffusion  along  a  2-D  grain  boundary 

(attached  as  Appendix  C) 


The  papers  published  so  far  in  this  and  preceding  ARO  contracts  arc  all  based  on  the  cluster  variation 
method  (CVM)'’  of  equilibrium  statistical  mechanics.  The  basic  approach  is  to  minimize  the  free  energy 
(actually  the  grand  potential)  with  respect  to  the  basic  variables  to  derive  the  equilibrium  state. 


The  CVM  has  been  extended  to  non-equilibrium  processes  by  treating  the  time  axis  as  the  fourth 
space  axis.  The  extension  has  been  called  the  path  probability  method  (I’I’M)4  It  was  used  successfully 
m  treating  atomic  diffusion  in  crystals  '  ,  and  in  the  kinetics  of  phase  transitions. 


lit  the  attached  Appendix  C.  vve  have  used  the  PPM  to  study  atomic  diffusion  along  a  grain  boundary. 
A  chemical  potential  gradient  Yp  is  imposed  along  the  g.b.  to  induce  atomic  flux.  Qualitatively  speaking, 
the  ratio  of  the  flux  and  Vp  gives  the  diffusion  coefficient.  1'he  treatment  is  based  on  the  equilibrium 
formulation  in  Appendix  B  and  works  on  the  impurity  g.b.  case. 


Taking  advantage  of  the  finding  in  the  following  section  4,  the  diffusion  is  formulated  in  the  skeletal 
crystalline  lattice  only  (ignoring  the  rest  of  the  DSC  lattice  points)  using  the  pair  as  the  basic  cluster.  This 
is  drastically  simpler  than  using  a  nine  point  QS  as  the  basic  cluster  of  the  diffusion  formulation. 


'Die  work  has  a  novel  technical  contribution  In  previous  PPM  formulations6,7  of  diffusion,  only  the 
bulk  diffusion  was  treated,  and  thus  the  gradients  of  state  variables  arc  small  under  the  imposed  Vp.  In 
the  present  case,  the  gradients  of  states  variables  arc  not  small  near  the  g.b.  It  is  shown  in  Appendix  C 
that  even  when  the  space  gradients  of  state  variables  in  equilibrium  are  large,  atomic  fluxes  can  be 
formulated  when  small  Vp's  arc  imposed. 


V  I  heory  of  diffusion  along  a  2-D  grain  boundary  (attached  as  Appendix  C) 


4 


When  a  reasonable  set  of  kinetic  parameters  arc  used,  numerical  computations  lead  to  the  results  that 
the  diffusion  is  many  orders  of  magnitude  faster  near  the  g.b.  than  inside  the  bulk.  This  result  holds  even 
when  the  impurities  are  adsorbed  to  or  repelled  from  the  boundary.  The  diffusion  coefficient  is  roughly 
proportional  to  the  square  of  the  local  vacancy  concentration. 

The  diffusion  coefficient  along  a  line  parallel  to  the  g.b.  follows  the  Arrhenius  relation  nicely,  and  the 
activation  energy  can  be  determined  for  each  line  at  a  certain  distance  from  the  g.b.  The  activation  energy 
is  smaller  near  the  g.b.,  corresponding  to  the  accumulation  of  vacancies  near  the  g  b.  center. 


.V  Theory  of  diffusion  along  a  2-0  grain  boundary  (attached  as  Appendix  (') 


4.  2-D  grain  boundary  between  ordered  phases 

(attached  as  Appendix  D) 


Intending  the  work  in  Appendix  B,  a  2-D  g.b.  between  two  ordered  phases  is  studied.  The  DSC 
lattice  of  Appendices  A  and  B  form  the  basis  of  tire  analysis,  and  the  stable  cry  stalline  state  is  the  square 
lattice  of  the  same  shape.  In  this  paper,  the  skeletal  square  lattice  is  ordered  when  the  system  is  close  to 
the  stoichiometric  AB  alloy. 


Two  approaches  were  used  in  the  formulation.  One  is  based  on  the  DSC  and  uses  the  quadruple 
square  (QS)  mentioned  in  Section  2  as  the  basic  cluster.  T  he  other  works  only  with  the  skeletal  square 
lattice  and  uses  a  pair  of  lattice  points  in  this  lattice  as  the  basic  cluster.  A  significant  lindiim  is  that  the 
two  approaches  give  identical  numerical  results  for  low  temperatures  including  the  entire  ordered  region 

l:or  the  coherent  boundary  of  the  present  study,  the  symmetric  boundary  and  the  asymmetric 
boundary  arc  distinguished.  I 'ruler  a  certain  choice  of  interaction  parameters,  the  stable  boundary  is 
symmetric  in  low  temperatures  but  switches  to  asymmetric  in  high  temperatures. 

I  he  bulk  shape  of  the  skeletal  crystalline  lattice  is  the  square  lattice.  However,  near  the  g.b.  the  shape 
is  distorted  and  thus  the  pair  cluster  is  insufficient  in  taking  into  account  the  frustration  effect.  It  was 
found  that  an  additional  triangle  cluster  near  the  coincidence  lattice  site  can  remedy  the  entropy  expression 
and  leads  to  the  correct  kln2  value  at  the  OK  limit. 


T  2-D  train  boundary  between  ordered  phases  (attached  as  Appendix  D) 


5.  Transition  layer  in  a  lattice-gas  model  of  a 
solid-melt  interface  (attached  as  Appendix  L) 
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Appendix  F 


Complete  summary  of  the  author's  work  on  boundary  structure 
related  to  the  present  contract. 


Since  the  present  final  report  closes  the  scries  of  projects  on  boundary  structure  supported  by  the 
IJ.S.  Army  Research  Office,  and  a  new  page  of  the  boundary  studies  is  about  to  open,  it  is  appropriate 
to  summarize  the  entire  work  done  by  the  author  in  the  subject  of  interest. 


1.  CVM  and  PPM 


These  two  methods  form  the  basis  of  the  boundary  studies. 


(K51)  R.  Kikuchi,  "A  theory  of  cooperative  phenomena,"  Phys.  Rev.  XI  (1951)  9SX.  This  is  the  first 
full  paper  on  the  CVM,  and  is  still  a  good  introduction. 


(KB67)  R.  Kikuchi  and  S.(i.  Brush,  "Improvement  of  the  cluster-variation  method,"  .1.  ('hem  Phys. 
47  (1967)  195.  Tor  a  2-D  systin,  1-D  long  clusters  arc  enough  to  obtain  good  accuracy. 


(K74)  R.  Kikuchi,  "Superposition  approximation  and  natural  iteration  calculation  in  cluster-variation 
method,"  .1.  Chcm.  Phys.  60  (1974)  1071.  The  NI  method  is  an  important  computational  device  to  use 
in  the  CVM. 


(K66)  R.  Kikuchi,  "The  path  probability  method,"  Prog.  Thcor.  Phys.,  Suppl.  No.  35,  p.  I  (1966) 
Still  a  good  introduction  to  the  PPM. 
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2.  Karly  papers 


(K6(V)  R.  Kikuchi,  "Statistical  dynamics  of  boundary  motion,"  Annals  of  I’hys.  11  (I960  328.  This 
is  the  first  paper  by  the  author  on  the  boundary  structure.  It  is  interesting  that  the  work  was  on  kinetics. 
I'hc  Ising  model  boundary  between  the  +  and  -  domains  is  moved  perpendicular  to  itself  under  an 
external  magnetic  field.  The  velocity  of  the  motion  is  derived  as  an  eigenvalue  of  a  differential  equation. 
The  formulation  was  done  using  the  point  approximation  of  the  PPM. 


((.'KM')  J.W.  Cahn  and  R.  Kikuchi.  "Theory  of  domain  walls  in  ordered  structures  I  Properties  at 
absolute  zero."  .1.  I’hys.  ('hem.  Solids  20  (1961)  94. 


(KC62)  R.  Kikuchi  and  J.W  Cahn,  'Theory  of  domain  walls  in  ordered  structures  II.  Pair  approxi¬ 
mation  for  nonzero  temperatures,  '  .1.  Plus,  (  hem  Solids  23  (1962)  137. 


(CK66)  JAV  Cahn  and  R.  Kikuchi,  Theory  of  domain  walls  in  ordered  structures  III  (  fleet  of 
sutstitutional  deviations  from  stoichiometry."  I  Plus  Chem.  Solids  27  (1966)  1305,  In  these  three  papers, 
the  CVM  was  used.  However,  since  the  M  method  in  (K74)  below  was  not  known,  only  the  pair 
approximation  was  applied  for  an  antiphase  boundary  (APB)  in  bee. 


.3.  Ilic  scalar  product  (SP)  formulation 


The  excess  free  energy  of  the  boundary  is  calculated  using  only  the  bulk  properties  of  the  two  sides. 


(K67)  R.  Kikuchi,  "Cooperative  phenomena  in  the  triangular  lattice,"  .1.  Chem.  I’hys.  47  (1967)  1664. 
This  paper  shows  that  the  “central  hump"  of  the  free  energy  curve  disappears  as  we  improve  approximation 
in  the  CVM  hierarchy.  This  conclusion  makes  us  question  the  validity  and  or  accuracy  of  the  CVM 
treatments  of  the  boundary  structure,  for  example  in  the  papers  in  2  above. 


((  \V7|)H  B  Clayton  and  (i.W  Woodbury,  Jr,  J.  Chem  I’hys.  5s  (1971)  (S95.  Ibis  paper  shows 
how  to  calculate  the  excess  free  energy  of  a  boundary  using  the  bulk  properties  only  of  the  two  sites  1  heir 
method  is  used  in  the  subsequent  papers. 
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(K721, 11,111)  R.  Kikuchi,  "Boundary  free  energy  in  the  lattice  model.  1.  General  formulation,"  .1. 
('hem.  Pins.  57  (1072)  777;  II.  Applications  of  the  general  formula,".!,  ('hem.  Phvs.  57  (1072)  782;  "III. 
Solution  of  the  paradox,"  .!.  ('hem.  Phys.  57  (1072)  787.  These  three  papers  present  detailed  discussions 
of  the  SP  method,  and  the  interpretation  of  the  CVM  treatment  in  2. 

(K72IV)  R.  Kikuchi,  "Phase  transition  within  a  phase  boundary,"  .1.  ('hem.  Phys.  57  (1072)  46.72. 
I 'sing  the  SP  method,  it  is  shown  that  an  order-disorder  type  phase  transition  occurs  within  a  phase 
boundary,  when  the  boundary  itself  is  2-D. 

(K7.7)  R.  Kikuchi,  "Phase  transition  within  a  phase  boundary  II,"  Scripta  Met  7  (1077)  275.  This 
compares  works  of  several  different  groups  on  the  subject. 

(K76)  R.  Kikuchi,  "Natural  iteration  method  and  boundary  free  energy,"  .1.  ('hem.  Phys.  65  (2076) 
4545.  I  rom  this  paper  on,  the  boundary  work  was  supported  by  the  I  S.  Army  Research  OITicc.  The 
paper  uses  the  NI  method  to  calculate  the  boundary  free  energy  formulated  by  the  SP  method. 

(K77a)  R.  Kikuchi,  The  scalar- product  experssion  of  boundary  free  energy  for  long-range  interacting 
systems,"  .1.  ('hem.  Phys.  68  (1077)  110.  The  SP  formulation  is  extended  to  the  case  when  the  interaction 
is  not  limited  to  the  nearest-neighbor  pairs. 

4.  Introduction  of  the  NIM 

The  natural  iteration  method  (NIM)  was  devised  in  (K74)  quoted  in  1  above.  The  NIM  made  it 
possible  to  calculate  boundary  structure  with  basic  clusters  larger  than  the  pair. 

(K77b)  R.  Kikuchi,  "Structure  of  phase  boundaries,"  .1.  ('hem.  Phys.  66  (1077)  2252.  In  addition  to 
the  NIM.  a  special  technique  adopted  from  the  work  of  Weeks  and  Gilmer  is  used  to  calculate  <  llf)> 
boundary  of  the  bcc  King  model  using  a  tetrahedron  as  the  basic  cluster.  Results  arc  compared  with  the 
SP  treatment  of  the  same  boundary. 

5.  I  ce  boundary  (ll’B  and  APB)  using  a  tetrahedron 
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(CK79)  J.W.  Cahti  and  R.  Kikuchi,  "Ground  slate  of  fee  alloys  with  multiatom  interactions,"  Aeta 
Met.  27  (1079)  1320. 


& 


(KC79)  R.  Kikuchi  and  J.W.  Cahn,  'Theory  of  interphase  and  antiphaes  boundaries  in  fee  alloys," 
Aeta  Met.  27  (1970)  1337.  These  two  papers  form  a  connected  pair.  II’R  and  APB  of  CujAu  arc 
calculated  using  the  tetrahedron  as  the  cluster.  The  bulk  phase  is  the  best  analytical  formulation  calculated 
by  van  Baal,  Golosov,  Kikuchi,  and  de  Fontaine  using  the  CVM. 
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(CKR5)  J.W.  (,'ahn  and  R.  Kikuchi,  "Transition  layer  in  a  lattice-gas  model  of  a  solid-melt  interface, 
I’hys.  Rev.  B31  (1985)  4300.  This  is  the  paper  in  Appendix  F,  of  this  report. 


6.  2-1)  grain  boundary  studies. 


This  is  the  last  category  of  the  boundary  studies  done  by  the  author  so  far. 
Appendices  A  through  D  of  this  report. 


The  papers  arc  in 
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